Kolmogorov Spectrum of Quantum Turbulence 
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There is a growing interest in the relation between classical turbulence and quantum turbulence. 
Classical turbulence arises from complicated dynamics of eddies in a classical fluid. In contrast, 
quantum turbulence consists of a tangle of stable topological defects called quantized vortices, and 
thus quantum turbulence provides a simpler prototype of turbulence than classical turbulence. In 
this paper, we investigate the dynamics and statistics of quantized vortices in quantum turbulence 
by numerically solving a modified Gross-Pitaevskii equation. First, to make decaying turbulence, 
we introduce a dissipation term that works only at scales below the healing length. Second, to 
obtain steady turbulence through the balance between injection and decay, we add energy injection 
at large scales. The energy spectrum is quantitatively consistent with the Kolmogorov law in both 
decaying and steady turbulence. Consequently, this is the first study that confirms the inertial range 
of quantum turbulence. 
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PACS numbers: 67.40.Vs, 47.37.+q, 67.40.Hf 

I. INTRODUCTION 

The study of turbulence has very long history, going 
back at least to Leonard da Vinci, and the turbulent 
state in fluid dynamics is now intensely studied in vari- 
ous fields including mathematics, and the experimental, 
theoretical, and computational branches of physics 0,0- 
Until recently, only classical turbulence (CT) has been 
studied. Classical incompressible Newtonian fluids are 
usually described using the Navier-Stokes equation 



0_ 

dt 



v(x, t) + v(x, t) ■ Wv(x, t) 
= — VP(as, t) + vW 2 v(x, t) 



V • v(x,t) = 0, 



(la) 
(lb) 



where v(x,t) is the velocity of fluid, P(x,t) is the ex- 
ternal pressure, p is the density of fluid and v is the 
kinematic viscosity. Time development of v(x, t) is dom- 
inated by eq. Ijlafl . while eq. (|lb|) confirms incompress- 
ibility of fluid and determines P(x, t). Flow of this fluid 
can be characterized by the ratio of the second term of 
the left-hand side of eq. l|la[) . hereafter the inertial term, 
to the second term of the right-hand side, hereafter the 
viscous term. This ratio is the Reynolds number 



v 



(2) 



Here v and L are the characteristic velocity of flow and 
the characteristic scale, respectively. When v increases 
to the point that the Reynolds number exceeds a critical 
value, the system changes from a laminar to a turbulent 
state in which the flow is highly complicated with many 
eddies. These eddies have a high concentration of vortic- 
ity u(x,t) = V x v(x,t). Owing to its complexity, the 
dynamics of the turbulent state is very difficult to inves- 
tigate. Nevertheless, G. I. Taylor, a pioneer of modern 
studies of turbulence, built a theoretical and experimen- 
tal framework for turbulence by introducing the idea in 



which turbulence is homogeneous and isotropic at smaller 
scales than the system size L Based on this idea, 
Kolmogorov proposed a self-similar, equilibrium state of 
fully developed turbulence together with a statistical law 
of turbulence called the Kolmogorov law 0, . The Kol- 
mogorov law was confirmed through many experimen- 
tal and theoretical studies, which led to an explosion in 
our understanding of CT. Today, numerical simulation 
of the Navier-Stokes equation is an indispensable tool 
for understanding CT because it is almost impossible to 
study the equation analytically. The rapid development 
of computers with the resulting improvements of large- 
scale simulations has accelerated our understanding of 
CT. However, two major problems have remained. One 
is that a large system size L and a sufficient spatial res- 
olution are required for analyzing well-developed CT at 
high Reynolds numbers R; a problem that needs higher- 
powered computers for its solution. The other problem is 
the identification of eddies in CT: there is no global way 
to identify classical eddies. There are several variables to 
identify eddies, such as the pressure, the velocity field, 
and the vorticity [||, but the identification of eddies in 
CT depends on the choice of these variables and their 
thresholds of identification. As a result, the dynamics of 
eddies, for example, the cascade process of eddies in CT, 
becomes only conceptual. Consequently, the existence of 
the cascade process, the self-similarity of the dynamics of 
eddies, and how they are related to the statistics of CT 
are not well understood. 

Recently, in the field of low temperature physics, a 
similarity has been found between CT and turbulence in 
superfluid 4 He. As a result, there has been much interest 
in studying superfluid turbulence (ST) as a way to better 
understand CT. 

The physics of ST in liquid 4 He is one of the most im- 
portant topics in low temperature physics . Liquid 4 He 
enters the superfluid state at 2.17 K. Below this temper- 
ature, the hydrodynamics is usually described using the 
two-fluid model in which the system consists of an invis- 
cid superfluid and a viscous normal fluid. Early experi- 
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mental studies on the subject focused on thermal coun- 
tcrflow in which the normal fluid and superfluid flow in 
opposite directions. Their flow is driven by the injected 
heat current, and it was found that the superflow be- 
comes dissipative when the relative velocity between two 
fluids exceeds a critical value [jj . Feynman proposed that 
this is a superfluid turbulent state consisting of a tangle of 
quantized vortices |8(- Vinen later confirmed Feynman's 
proposition experimentally by showing that the dissipa- 
tion comes from the mutual friction between vortices and 
the normal flow 0]. After that, many experimental stud- 
ies were done on ST in thermal counterflow systems . 
In particular, Schwarz clarified the picture of ST having 
tangled vortices by doing a large-scale numerical simu- 
lation of the quantized vortex-filament model in thermal 
counterflow However, as thermal counterflow has 

no analogy with conventional fluid dynamics, it has not 
helped us understood the relation between ST and CT. 

More recently, experiments have been done on ST that 
do not involve thermal counterflow. Such studies opened 
a new stage in studies of ST. Maurer et al. studied a 
ST that was driven by two counter-rotating disks in 
superfluid 4 He at temperatures above 1.4 K. They mea- 
sured the pressure and converted it to the local velocity 
fluctuation, and found that the energy spectrum was con- 
sistent with the Kolmogorov law. Stalp et al. measured 
the attenuation of second sound [13j in superfluid 4 He at 
temperatures above 1 K, and found that inferred decay 
of grid turbulence was consistent with the Kolmogorov 
law. These experiments are consistent in the sense that 
they show similarities between ST and CT. This similar- 
ity can be understood using the idea that the superfluid 
and the normal fluid are likely to be coupled together 
by the mutual friction between them and thus to behave 
like a conventional fluid [l4|. However, the normal fluid 
is negligible at very low temperatures, and this brings up 
an important question: even without the normal fluid, is 
ST still similar to CT? Such a fluid is a pure quantum 
fluid and its turbulence can be called quantum turbulence, 
hereafter QT. 

To address this question, we consider the statistical 
law of CT 2]. For fully developed turbulence of an in- 
compressible classical fluid, the energy is injected into 
the fluid at scales comparable to the system size L in the 
energy-containing range, and in the inertial range, this 
energy is transferred to smaller scales without being dis- 
sipated. In this range, the system is locally homogeneous 
and isotropic, which leads to the statistics of the energy 
spectrum known as the Kolmogorov law 



E(k) = Ce'^k-^ 3 . 



(3) 



Here the energy spectrum E(k) is defined as E = 
J dk E{k), where E is the kinetic energy per unit mass 
and k is the wave number from the Fourier transfor- 
mation of the velocity field. The energy transferred to 
smaller scales in the energy-dissipative range is dissipated 
through the viscosity of the fluid with the dissipation rate 
e in eq. which is equal to the energy flux n in the 



inertial range. The Kolmogorov constant C is a dimen- 
sionless parameter of order unity . 

The inertial range is thought to be sustained by the 
self-similar Richardson cascade in which large eddies are 
broken up into smaller ones with many reconnections. In 
CT, however, the Richardson cascade is not completely 
understood because it is impossible to definitely identify 
each eddy - the eddies keep nucleating and annihilating 
due to the viscosity. In contrast, quantized vortices in 
a quantum fluid are definite and stable topological de- 
fects. At low temperatures, a macroscopic number of 
particles in a Bose system, such as liquid He, occupy a 
single particle ground state. In this case, matter waves of 
Bose particles form a coherent scalar macroscopic wave 
function &(x, t) = f(x, t) exp[i0(;c, t)], which is known as 
Bose-Einstein condensation. The dynamics of the macro- 
scopic wave function is governed by the Gross-Pitaevskii 
(GP) equation [H [13 

i t) = [-V 2 -n + g\$(x, t)\ 2 Mx, t). (4) 

The CP equation is obtained from the Hamilton- Jacobi 
equation for the Hamiltonian 



H = I dx 



\^(x,t)\ 2 - ^(x,t)\ 2 + ^(x,tt], (5) 



where fj, is the chemical potential, g is the coupling con- 
stant, |<E>(;c,i)| 2 = f(x,t) 2 is the condensate density, 
and the superfluid velocity v(x,i) is given by v(x,t) = 
2V<j>(x,t). The vorticity u>(x,t) — V x v(x,t) van- 
ishes everywhere in a single-connected region of the fluid 
and thus, all rotational flow is carried by quantized vor- 
tices. In the core of each vortex, $(x,t) vanishes so 
that the circulation § v ■ ds around the core is quan- 
tized by 47r. The vortex core size is given by the healing 
length £ = l/f^/g. Quantized vortices at finite temper- 
atures can decay through mutual friction with the nor- 
mal fluid at any scale. At very low temperatures the 
vortices can decay by 1) the emission of compressible 
excitation through vortex reconnections and 2) a 

process in which vortices reduce through the Richardson 
cascade process and eventually change to compressible 
excitations. These are the only two decay mechanisms 
at very low temperature. In any case, dissipation occurs 
only at scales below the healing length; for larger scales, 
we can obtain the turbulent state at high Reynolds num- 
bers that is free from dissipation by quantized vortices. 
Therefore, QT at very low temperatures gives the high 
Reynolds number and the real Richardson cascade of def- 
inite quantized vortices, and thus can become an ideal 
prototype to study statistics such as the Kolmogorov law 
and the Richardson cascade in the inertial range of CT. 

There are two formulations for studying the dynamics 
of quantized vortices in QT. One is the vortex- filament 
model [TT| . and the other is the CP model. By using 
the vortex-filament model with no normal fluid compo- 
nent, Araki et al. studied numerically a vortex tangle 
arising from Taylor-Green flow, and obtained an energy 
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spectrum consistent with the Kolmogorov law |19|. By 
eliminating the smallest vortices, which had a size compa- 
rable to the numerical space resolution, they introduced 
dissipation into the system and realized decaying turbu- 
lence. Nore et al. used the GP equation to numerically 
study the energy spectrum of QT [2£| • They showed that 
the kinetic energy consists of a compressible part due to 
compressible excitations and an incompressible part com- 
ing from quantized vortices. Compressible excitations of 
wavelength smaller than the healing length are created 
through vortex reconnections or through the disappear- 
ance of small vortex loops flU l2ll| , so that the incom- 
pressible kinetic energy changes to compressible kinetic 
energy while conserving the total energy. The spectrum 
of the incompressible kinetic energy is consistent with 
the Kolmogorov law in a short period |22(. However, 
the consistency is broken in the late stage when many 
compressible excitations are created [lg, |2l|, |23|, |24| and 
hinder the Richardson cascade process, even on the large 
scales. 

Our approach here is to introduce a dissipation term 
that works only on scales smaller than the healing length 
£. This dissipation term can account for dissipation pro- 
cesses such as the scattering between the particle reser- 
voir and the emitted compressible excitations that can 
occur only at scales smaller than £ even at T = K 
[Til l2ll | . This dissipation removes only short- wavelength 
excitations (not vortices), thus preventing the vortices 
from receiving the excitation energy. Compared to the 
usual GP model, this approach enables us to more clearly 
study the relation between QT and CT. The Richardson 
cascade process at the large scales must be independent 
of the details of the dissipative mechanism at the small 
scales. 

We use our modified GP equation to investigate two 
types of QT. For decaying turbulence, we start from ran- 
dom configurations of the phase of the wave function and 
study the decay. We provide here a more detailed study 
of this type than that discussed in our previous paper 
25]. For steady turbulence (the other type), we assume 
that energy is injected into the system at a large scale by 
a moving random potential. For both types, compressible 
excitations are highly suppressed compared to quantized 
vortices, and the spectrum of the incompressible kinetic 
energy is found to be consistent with the Kolmogorov 
law, which verifies that QT is similar to CT. We also 
show here that steady turbulence with the moving ran- 
dom potential clarifies the energy-containing range, the 
inertial range and the energy-dissipative range. By calcu- 
lating the energy dissipation rate and the energy flux in 
the wave number space, we confirmed the following sce- 
nario. First, vortices are nucleated by the moving ran- 
dom potential in the energy-containing range. Second, 
the vortices decay into smaller vortices in the inertial 
range through the Richardson cascade process. Finally, 
the smallest vortices change to compressible excitations, 
which are then removed by the introduced dissipation. 
Although this scenario of quantized vortices in QT is very 



similar to that believed about eddies in CT, it is much 
clearer, because quantized vortices are definite topologi- 
cal defects. 

This paper is organized as follows. In Sec. [HJ we re- 
view our model of the modified GP equation that includes 
small-scale dissipation. After discussing the numerical 
procedure of our model in Sec. IIHI we show numerical 
simulations of i) a straight vortex line in Sec. IIVI ii) de- 
caying turbulence in Sec. and iii) steady turbulence 
in Sec. I VII Section IV 111 is devoted to our conclusions. 

II. GP EQUATION WITH SMALL-SCALE 
DISSIPATION 

In this section, we discuss the dynamics of the origi- 
nal GP equation Q with quantized vortices and why the 
compressible short-wavelength excitations make it diffi- 
cult to study QT. Then we propose a modified GP equa- 
tion with small-scale dissipation. The modified GP equa- 
tion gives a clear criterion for the inertial range of QT 
and has a form similar to that of the Navier-Stokes equa- 
tion with large effective Reynolds number. 

The dynamics of quantized vortices was analyzed by 
Koplik et al [2(|. and Ogawa et al. [2l| by using the 
GP equation Through the Madelung transforma- 

tion &(x,t) = f(x, t) exp[icj)(x, t)}, the GP equation (@J 
becomes 

f{x, t)^(x, t) =V 2 f{x, t) - f(x, t)V0(x, tf 

-[gf{x,tf-Af{x,t). (6) 

At scales larger than £, we can neglect the curvature of 
the amplitude of the wave function near quantized vor- 
tices Vf{x,t). With this assumption, the equation for 
the vorticity uj(x,t) = V x v(x,t) — 2V x X7(f>{x,t) be- 
comes 

^{x,t) + ^V xVv{x,tf = Q. (7) 

Thus, the dynamics of quantized vortices depends only 
on the fluid velocity v(x,t), which is similar to that of 
isolated eddies in a perfect classical fluid. Except for the 
dynamics of vortices described in eq. 10, there are two 
non trivial processes of quantized vortices; one is the re- 
connection of two vortices and the other is the disappear- 
ance of small vortex loops. Figure shows an example 
of vortex reconnection given by the numerical calcula- 
tion of eq. starting from two straight vortex lines in 
a skewed orientation. These two vortices approach and 
become locally anti-parallel, and then the reconnection 
occurs. After the reconnection, the vortices tend to be 
bent, which allows smaller structure to form, thus sus- 
taining the Richardson cascade process. Ogawa et al. 
show that this process does not violate Kelvin's circula- 
tion theorem; moreover, the process becomes irreversible 
due to the emission of compressible excitations that have 
a wavelength that is smaller than £. In the scales, its 
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(a) (b) (c) (d) 

FIG. 1: The reconnection of two quantized vortices starting from a skew position of two straight vortex lines, (a) Initial state, 
(b) Just before the two vortices connect, (c) Connection of two vortex lines, (d) After the vortex lines separate in a newly 
connected arrangement, also called reconnection. Contours in all figures are high concentrations of vorticity (98% of maximum 
vorticity) . 



dynamics is strongly disturbed by compressible short- 
wavelength excitations described by Vf(x,t), which is 
neglected in eq. Q. In turbulent states, such excitations 
are rapidly created through reconnections and disappear- 
ances of small vortex loops. Therefore, the dynamics 
of QT is affected by many compressible excitations and 
one cannot see the proper dynamics of quantized vortices 
from eq. Q by itself. 

To obtain QT free from compressible short-wavelength 
excitations, we have to introduce a dissipation term in 
the GP equation. First, we discuss a possible dissipation 
mechanism in the GP equation. Next, the dissipation 
term is arranged so that it can remove short-wavelength 
excitations. We assume that the system is described by 
&(x,t) such that the energy and particles are exchanged 
with a particle reservoir. The particle reservoir is a ther- 
modynamic environment that lets the chemical potential 
of the system equal that of the reservoir. The interaction 
with the particle reservoir is provided by an imaginary 
term in the Hamiltonian eq. @ 



H = jdx |V$(a;,t)| 2 - ^(x,t)\ 2 

+ ^(x,t)\ i -ir(x,t)\<s>(x,t)\ 2 

which yields the dissipative GP equation 



(8) 



= [-V 2 +g\^ tl (x,t)\ 2 ~iT(x,t)}^ ll (x,t). (9) 



Here the wave function <E> M (a:,£) includes the chemical 
potential /i through the gauge transformation ^(a;, t) = 
§(x, t) exp(— i/i£). Since T(x,t) comes from the differ- 
ence of the chemical potential, we can write T(x, t) — 
j(x,t)(fj, — /i pr ) with the chemical potential of the par- 
ticle reservoir being /i pr . We assume that the system is 
nearly in equilibrium with the particle reservoir and that 
T(x,t) is proportional to the difference in the chemical 



potentials. Using the approximation 
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(10) 



and ji ~ /i pr , eq. (JSJ) becomes 



d 

[i-70M)]^$„0M) 



= [-V 2 +s|* P (*,t)| 2 +iM7(*,0]*M( a! >t). (11) 

Substituting <5> p (a?,t) = $(x, t) exp(— ifit) into eq. JTTf l . 
we finally obtain the modified GP equation with the ef- 
fective dissipation term j(x,t) p7| 

[i-7(M)]J^0M) 

= {-V 2 -[i(t)+g\<S>(x,t)\ 2 ]<S>(x,t). (12) 

Introduction of j(x, t) conserves neither the energy nor 
the number of particles. However, for studying the hy- 
drodynamics of turbulence, it is realistic to assume that 
the number of particles is conserved. Hence, we consider 
the time dependence of the chemical potential so 
that the total number N — j dx \<fr(x,t)\ 2 can be con- 
served. The approximation of eq. l|l(Jf) gives the following 
form for the chemical potential 

(i(t) = ^ J dz<r(^i)[-V 2 +g|$(^t)| 2 ]<I>(:r,i). (13) 

We consider the spatial and temporal dependence 
of j(x, t) so that this term can dissipate the short- 
wavelength excitations. The Fourier transformation of 
the GP equation i|12|) can be written as 



[i-7(fc)]JU(M) 

= [fc 2 - M (t)]l>(M) 



L 6 



x $(fc - fei + k 2l t). 



(14) 
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Here &(k,t) is the spatial Fourier component of <f>(a;,i) 
and L is the system size. 7(fe) is the wave number 
representation of the dissipation term. For the dis- 
sipation term to remove only compressible excitations 
of wavelength smaller than £, it should have the form 
7(fe) = r yoO(k — 2ir/£) with the step function 9. This 
form is based on the idea that only compressible short- 
wavelength excitations interact with the particle reser- 
voir. Because system is dissipationless at scales exceed- 
ing £, we can investigate proper vortex dynamics at large 
scales without considering dissipation. 

Next we consider the Reynolds number of QT. 
By using the Madelung transformation Q(x,t) = 
f(x, t) exp[i(f>(x, t)], eq. |Q has the hydrodynamic form 
in the first order of 7(0;, t) 

d 

f{x,t)—6(x,t) 

= j(x,t)[2Vf(x,t) ■ V0(x,t) + f(x,t)V 2 <j>(x,t)] 
+ V 2 f(x,t)-f(x,t)[V<b(x,t)] 2 
-[gf(x,t) 2 ^fi(t)]f(x,t). (15) 

By neglecting the term Vf(x,t) and Vj(x,i) (because 
they are ineffective at scales exceeding £), we obtain 
the Navier-Stokes-like equation of the superfluid veloc- 
ity v(x,t) = 2V<p(x,t), 



^-v(x,t) + ^Vv(x,t) 2 

= 7(aj, t)V[V • v(x, t)] - 2V[s/(£c, tf - /z(i)]. (16) 



A comparison with the Navier-Stokes equation llajl 
shows that the second term on the left-hand side and the 
first term on the right-hand side of eq. Ijltjfl can be re- 
garded as the inertial and dissipation terms, respectively, 
of the GP equation, and also that -f(x, t) acts as the vis- 
cosity in a conventional fluid. Following this comparison, 
the Reynolds number Rqt of QT is 



Rqt 



vL 

— 1 

7 



(17) 



where v is the mean velocity of QT. Here 7 is the av- 
erage of 7(a3,t) over space and time. The average over 
the whole volume corresponds to the zero-wave number 
component of the Fourier transformation, namely 7 ~ 
7(fe = 0). However, the average over the sub- volume of 
dimension smaller than £ becomes 7 ~ 7(|fe| > 27r/£) = 
7o- Consequently, we obtain a scale-dependent Reynolds 
number Rqt which is basically different from that of clas- 
sical turbulence. At scales exceeding £, the Reynolds 
number Rqt thus becomes infinite, which means that 
scales with k < 27r/£ can be regarded as the inertial 
range of QT. At scales below £, Rqt has a finite value 
vi/70. 



III. NUMERICAL PROCEDURE 

To solve the GP equation numerically with high accu- 
racy, we use the pseudo-spectral method in space with 
periodic boundary conditions in a box with spatial res- 
olution containing A gr id = 256 3 grid points. In this 
box, $(x,i) is discretized to 3>[lAx, mAt], $(fc,i) is dis- 
cretized to $[ZAfc, mAt], /j,(t) is discretized to fi[mAt], 
and 7(fe) is discretized to 7[ZAfc]. The first two are re- 
lated through the Fourier transformation 

$\lAx, mAt] = V $[i'Afc, mAt] 

N„ id Ax A ^ 



V 



x exp[i(Z' • l)AkAx], 



(18) 



where integers m and Z = [l x ,ly>lz] satisfy m > and 
< l x ,lyJz < 256. Also, the resolutions Aar, Afc = 
2n / L and At are for space, wave number, and time. To 
obtain i>[ZAfc, (m + I) At] from $[ZAfc,mAt], we have to 
calculate the right-hand side of eq. I|14|). However, it 
takes a long time (0(AT| rid )) to calculate the second term, 
so we instead calculate h(x,t) = g\$(x, t)\ 2 $(x, t) and 
its Fourier transformation h(k, t). Then the GP equation 
becomes 



[i- 7 (fc)]||.(fc,i) = [fc 2 - M (i)]$(fe,t) 



h(k,t). (19) 



Because h(x,t) is the functional of $>(x,t), we should 
calculate it again from $[ZAfc, (m + l)At] by using the 
Fast-Fourier- Transformation (FFT) [23] and the inverse 
FFT, for the next step $[ZAfc,(m + 2) At]. This pro- 
cedure reduces the computation time from 0(N^ lid ) to 
0(iV grid log iVgrid) of the FFT. At each time step, the 
wave function is normalized, so that the total number of 
particles N can be conserved, and the chemical potential 
is calculated using eq. ifHUl . 

In this paper, we consider the case of g = 1. For numer- 
ical parameters, we use a spatial resolution of Ax = 0.125 
and L = 32, where the length scale is normalized by 
the healing length £. The wave number resolution is 
Afc = 2ir/L ~ 0.196. The numerical time evolution of 
eq. <|19|) was calculated using the Runge-Kutta-Verner 
method [2^] with the time resolution At = 1 x 10~ 4 . 
The flowchart of the calculation is as follows: 

(i) Calculate h[lAx, mAt] from &[lAx,mAt]. 

(ii) Calculate h[lAk, mAt] by the FFT. 

(iii) Calculate <t[lAk,(m+ I) At] from eq. (jT9|) . 

(iv) Calculate <j>[ZAx, (m + l)At] using the inverse FFT. 

(v) <i>[ZAa;, (m + l)At] is normalized to conserve N. 

(vi) Calculate fi[(m + l)At] from eq. JT3J). 
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In the framework of our numerical calculation, we 
cannot deal with the large wave number region k > 
fe c = (27r/Ax,27r/Aa;,27r/Ax). Throughout the FFT, 
the component of the wave function $(fe > k c ,t) with 
such a large wave number cycles back to that with a small 
wave number, a numerical process that causes an alias- 
ing error. To reduce the aliasing error, we placed a large 
number of grid points A/grid in a box and set the value 
of k c such that the wave function is negligible near k c . 
We found that = 256 3 is large enough to neglect 

the aliasing error for our QT. We discuss this problem 
further in section Ivl 



IV. DYNAMICS OF A SINGLE STRAIGHT 
VORTEX LINE WITH DISSIPATION 

Before running the QT simulations, we investigated 
the effect of dissipation on the dynamics of a single vor- 
tex line. Here, we assume a flow, and then consider the 
resulting drag force due to dissipation. We use the scale- 
dependent dissipation term 7(fe) = r )o9(k — 2s7r/£) by 
introducing the parameter s, and study how the vortex 
dynamics depend on s. When s > 1 and the dissipation 
works only on scales below £, a vortex line moves free 
from these forces. Therefore, the dissipation term with 
s = 1 acts only on the short-wavelength excitations and 
does not affect the vortex dynamics at scales exceeding 
£. Thus s = 1 is suitable for this study. 

The steady solution of the GP equation with a single 
vortex line along the z axis is represented by 



$(x) = /(r)exp(kp), 



(20) 



where r = \/ x 2 + y 2 and ip = tan _1 (y/x) are the radius 
and angle in cylindrical coordinates. The amplitude of 
the wave function f(r) follows the equation 



d 2 /(r) | ld/(r) f(r) 
dr 2 r dr r 2 



+ f(r) - f(r) = 0, (21) 



and is obtained numerically by the Shooting-Runge- 
Kutta-Verner method |28| with the boundary conditions 
/(0) = and f(r — » oo) = 1. Starting from the single 
vortex solution (|20|l . we calculate the GP equation under 
a flow v c as 

= [-V 2 - /i + iv e ■ V + </|$(x, t)\ 2 }$(x, t). (22) 
The Madelung transformation of eq. lll'l'l) yields 
df(x,t) 



dt 



--v e • V/(aj, t) + -fix, t)f(x, t) v e • V0(as, t) 



j(x, t) 2 v e ■ Vf(x, t), 



(23) 



which is second order in j(x,t). The first term on the 
right-hand side is the inertial term for v e . The second 



and third terms of the right-hand side describe the mo- 
tion of the vortex line by the dissipation term ^{x, t). 
The direction of vortex movement is indicated by the 
sign of df(x,t)/dt in eq. (|23l) . as shown in Fig. [21 ( a )- 
As shown in this figure, the vortex moves from the side of 
increasing f(x,t) with time to that of decreasing f(x,t), 
that is, the positive to negative sides of df(x, t)/dt. Fig- 
ures El (b) and (c) show the directions of vortex motion 
caused by the second and third terms of the right hand 
side in eq. I)23|) respectively. The thin arrows in these 
figures show the direction of v e , thick arrows show the 
direction of V^>(cc, t) and — V/(a;, t) in Fig. 01(b) and (c), 
respectively, and the dashed lines show the boundary be- 
tween the positive and negative regions of v e ■ Vtfi(x,t) 
and —v c -\7f(x 7 1). For the case of Fig. [3(b), v e -V4>(x, t) 
becomes positive in upper side (and negative in the lower 
side), making the vortex move to the lower side. Conse- 
quently the second term in eq. I)23|) causes a drag force to 
act in the direction perpendicular to v c . By using similar 
considerations applied to v c ■ V<fi(x,t), we can conclude 
that — v c ■ Wf(x, t) makes the vortex move from right to 
left in Fig. 121(c), and the third term in eq. I|23|l describes 
the drag force parallel to v c . These two forces enable us 
to investigate the effect of dissipation j(x, t). 

We numerically solved the GP equation (|22[) under the 
flow v e = 0.5 with a dissipation of 70 = 1 and no dissipa- 
tion 70 = 0. The effect of dissipation can be determined 
by comparing the position of the vortex in two simula- 
tions. Here we let the dissipation have a scale-dependence 
as 7(fc) = 7o#(fc — 2s7r/£) with the parameter s, and mon- 
itor how the vortex shift depends on s. The position of a 
vortex can be identified by the singularity around which 
the phase rotates by 2n. Although the vortex just moves 
along v e for 70 = 0, the dissipation 70 = 1 changes its 
trajectory. The deviation Ax v (t) is divided into the part 

Ax^-{t) perpendicular to v c and the part Axv parallel 
to v c . Figure shows how the deviation depends on s. 
Both Aa;^- and Axl are negligible for s > 0.3 and they 
rapidly increase for s < 0.3. From these considerations, 
we conclude that the dissipation with s = 1 is adequate 
for the present study. 

At finite temperatures, vortices in QT dissipate at var- 
ious scales, which may be considered as the mutual fric- 
tion of superfluid 4 He. Also, the dependence of j(x, t) on 
s and 70 may correspond to the temperature-dependence 
of the mutual friction, but this possibility will not be 
determined here. 



V. DECAYING TURBULENCE 

In this section, we discuss a simulation of QT starting 
from a random phase profile |25j. With the dissipation 
term, the QT decayed gradually and fell into a quasi- 
steady turbulent state in which the vortices at large scales 
were nearly isotropic and uniformly. After showing that 
the decaying mechanism of QT is caused by the dissi- 
pation of compressible excitations, we obtain the energy 
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FIG. 2: (a) The amplitude of the wave function f(x,t) at t (solid line) and t + At (dashed line) across the line A-B in figure 

(b) and (c). (b) Movement of the vortex line in a cross section of the x — y plane, caused by the term v c ■ Vcj>(x, t) in eq. (IL'MI . 

(c) Same as (b) except for the term — v c ■ X7f(x,t). The dot at the center is the vortex core, the thick arrows are V(j)(x,t) (b) 
and — V/(x, t) (c), thin arrows are v e , and A and B represent the points in the positive and negative regions of v c ■ X7(f>(x,t) 
(b) and -v e ■ Vf(x,t) (c). 




(a) (b) 

FIG. 3: Components of the vortex deviation Ax v for a range 

of s values at t = 5. In the calculation, 70 = 1 and 70 = FIG. 4: Procedure for making a random phase, (a) Example 
under the velocity field v c = 0.5. of 4>o(x) with A = 4. (b) Example of an initial random phase 

in a cross-section <f>o(x,y,0). 



spectrum of such a decaying QT and show that it is con- 
sistent with the Kolmogorov law. 

To obtain a turbulent state, we start from an initial 
configuration in which the density f§ = 1 is uniform 
and the phase 4>o{x) has a random spatial distribution. 
Here, the random phase 4>o(x) is made by placing ran- 
dom numbers between —tt to tt at every distance A = 4 
and connecting them smoothly (Fig. by using a three- 
dimensional spline interpolation 28] . The initial veloc- 
ity v(x,t = 0) = 2W(f)o(x) given by the initial random 
phase is random; as a result, the initial wave function is 
dynamically unstable and soon produces fully developed 
turbulence with many quantized vortex loops. 

Before investigating the dynamics and statistics of QT, 
we discuss the numerical accuracy and the aliasing error. 
Because the aliasing error is maximized for 70 = 0, where 
the wave function with the high wave numbers is never 
dissipated, it would be enough to confirm that the alias- 
ing error is very small for a simulation with 70 = 0. We 
can estimate the aliasing error by calculating the dynam- 
ics with a larger wave number space, which decreases the 
aliasing error, and then comparing the result with the 
dynamics of the original space. Such a large wave num- 



ber space is realized by increasing the spatial resolution 
iVgrid = 512 3 (Ax = 0.0625). For the two spaces, we com- 
pare the total energy E(i), the interaction energy Ei nt (t), 
the quantum energy E q (t), and the kinetic energy E^i n {t) 



E(t) = 




r dx$*(x,t) -V 2 + 




E int {t) = 


9 

2N 


j dxf(x,t) 4 , 




E q (t) = 




r dx [\7f(x,t)} 2 , 




Ekin(t) = 




C dx [f(x,t)V(f>(x,t)} 2 


(24) 



In addition, we investigate the accuracy of the Runge- 
Kutta-Verner method for the time evolution by calcu- 
lating the dynamics with higher time resolution At = 
2 x 10~ 5 . For the comparison, we estimate the rela- 
tive errors F l2 {A) = \{{A) l - {A) 2 )/{A) l \ and F 13 (A) = 
\({A)i - (A) 3 )/(A) 1 \ at t = 14 for A = E(t), E int (t), 
E q (t) and -Ekin(i)- Here ()i means that the resolution 
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was (i) At = 1 x 10" 4 and 256 3 grid points, O2 means 
a resolution of (ii) At = 2 x 10" 5 and 256 3 grid points, 
and O3 means (iii) At = 1 x 10" 4 and 512 3 grid points. 
The comparison in Table [I] below shows that the rela- 
tive errors are extremely small, which confirms the small 
aliasing error and the high accuracy of the time evolu- 
tion, and allows us to use the original resolution of (i) 
above. The total energy is conserved to an accuracy of 

TABLE I: Dependence of F 12 and F 13 on E(t), E iat (t), £ q (t) 
and E hla (t). 





F 12 (A) 




Fi 3 (A) 




A -- 


= E(t) 2.5 x 10" 


12 


6.5 x 10" 


HI 


A = 


Bint(i) 3.9 x 10" 


12 


9.0 x 10" 


10 


A = 


= E q (t) 2.6 x 10" 


12 


7.1 x 10" 


10 


A = 


Eun(t) 5.2 x 10" 


12 


9.7 x 10" 


10 




(a) 




(b) 



FIG. 5: Time development of E(t), £ k m(£), E£ in (t), and 
^kin(^)- ( a ) Case with 70 = 0. (b) Case with 70 = 1 
(Kobayashi and Tsubota, Phys. Rev. Lett. 94, 065302-3, 
2005, reproduced with permission. Copyright (2005) by the 
American Physical Society). 



10" 10 . 

Next, we ran a simulation of QT with the dissipation 
70 = 1 to remove compressible short-wavelength excita- 
tions. The effect of this dissipation is shown in Fig. [5] 
We divide up the kinetic energy E^it) into the com- 
pressible part ££ in 0) = Jdx [{f(x,t)Vcj)(x,t)} c } 2 /N, 
which is due to compressible excitations, and the incom- 
pressible part E^Jt) = Jdx [{f(x,t)V(t)(x,t)y} 2 /N, 
which is due to vortices. Here {• • • } c denotes the com- 
pressible part V x {■ • ■ } c = and {■••}' denotes the 
incompressible part V-{- • • }' — 0. Numerically, the com- 
pressible part A c (x) and the incompressible part A 1 (x) 
of an arbitrary vector held A(x) are given by 

A c (x)=^ fc '^ (fc) fcexp(ifc-x) 
k 

A i (x) = A{x)- A c (x), (25) 

where A(k) is the Fourier component of A{x). Figure 
shows the time development of E(t), E-^i n (t), E^ in (t) 
and El in (t) for the cases of 70 = (a) and 70 = 1 
(b). Without dissipation, the compressible kinetic en- 
ergy E^ in (t) is increased in spite of the total energy con- 
servation (Fig. El (a)), a result that is consistent with the 
simulation by Nore et al. po| The dissipation suppresses 
E^ in (t) and thus causes E^ in (t) to be dominant in Fig. 
(b). This dissipation term acts as a viscosity term only 
at the small scales. At other scales, which are equivalent 
to the inertial range Afc < fc < 2n/£, the energy does not 
dissipate. 

To show that the result is consistent with the Kol- 
mogorov law, we calculate the spectrum of the incom- 
pressible kinetic energy E^ in (k,t) defined as E^ in (t) — 
J dk El in (k,t). Initially, the spectrum E^ in (k,t) signif- 
icantly deviates from the Kolmogorov power-law; how- 
ever, the spectrum approaches a power-law as the turbu- 
lence develops. We assume that the spectrum EL n (k, t) is 
proportional to fc"'' in the inertial range Afc < fc < 2-7r/£ 



and then determine the exponent r\ by making a straight 
fit to a log-log plot of the spectrum E^ a (k, t) only in the 
inertial range Afc < fc < 27r/£ at each time. The result- 
ing time development of rj is shown in Fig. [H] (a). This 
figure shows that the QT satisfies the Kolmogorov law 
El in (k,t) oc fc" 5 / 3 for time 4 < t < 10, the times when 
the system may be fully developed turbulence. We also 
calculate the energy dissipation rate e = —dEl in (t)/dt 
and compare the results quantitatively with the Kol- 
mogorov law. FigureHJ(c) shows that e is almost constant 
in the period 4 < t < 10, giving rj = 5/3, which means 
that the system can be considered to be a quasi-steady 
turbulent state with a constant Kolmogorov spectrum of 
e 2 / 3 fc"". 

Without the dissipation, that is, for 70 = 0, the QT 
did not satisfy the Kolmogorov law. This means that 
dissipating short-wavelength excitations are essential for 
producing the Kolmogorov law. Figure H3 (b) shows the 
time development of the exponent r\ for 70 = 0. In this 
case, r\ is consistent with the Kolmogorov law only in 
the short period 4 < t < 7. For t > 7, the dynamics 
of quantized vortices are affected by active compressible 
excitations. Such an effect of compressible excitations is 
also clearly shown in Fig. (d) . Compared with Fig. 
(c), e is unsteady and smaller than that of 70 = 1 and 
even becomes negative at times. When e is negative, the 
energy flows backward from compressible excitations to 
vortices, which may prevent the energy spectrum from 
satisfying the Kolmogorov law for t > 7. This result is 
consistent with the simulation by Nore et. al. |2£j Ac- 
cording to Fig. El (a), the agreement between the energy 
spectrum and the Kolmogorov law becomes weak at the 
later stage of t > 10, which may be due to the following 
reasons. In the period 4 < t < 10, the energy dissi- 
pation rate is constant and the energy spectrum agrees 
with the Kolmogorov law, which may support the find- 
ing that the system is in a quasi-steady turbulent state, 
and the Richardson cascade process works well at large 
scales. The dissipation is caused mainly by the removal 
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(a) (b) (c) (d) 

FIG. 6: Time dependence of the exponent of the spectrum r\ and the dissipation rate e. (a) rj for 70 = 1. (b) r\ for 70 = 0. 
(c) e for 70 = 1. (d) e for 70 = 0. In (a) and (b), the error bars are the standard deviation of the data from the fit. The line 
77 = 5/3 is shown to compare the results with the Kolmogorov law (Kobayashi and Tsubota, Phys. Rev. Lett. 94, 065302-3, 
2005, reproduced with permission. Copyright (2005) by the American Physical Society). 



of short-wavelength excitations that had been emitted at 
vortex reconnections. However, the system at the late 
stage t > 10, being no longer quasi-steady turbulent in 
large scales, has only small vortices after the Richard- 
son cascade process. Such a state is not regarded as a 
turbulent state, and thus the energy spectrums inconsis- 
tent with the Kolmogorov law of 77 = 5/3. Emissions of 
excitations through vortex reconnections hardly occurs, 
which greatly reduces the energy dissipation rate s. We 
can avoid such a disagreement by adding energy injection 
to the simulation. The result is steady turbulence, which 
is discussed in the next section. 

To confirm that the system is a fully developed turbu- 
lent state, we try to visualize the configuration of quan- 
tized vortices, the phase of which changes by 2tt around 
the core. Because this is different from vortices in the 
vortex-filament model [III fli^ , identification of vortices 
in the GP model is not straightforward. Nevertheless, 
the identification can be done using the following two 
methods. The first is based on the fact that vortex 
cores have nearly infinite vorticity; in this method, we 
determine the points with a huge vorticity \uj(x, t)\. We 
used this method in our previous study [25j- The other 
method is to directly search for the points around which 
the phase rotates by 2tt. For simplicity, we consider a 
straight vortex line along the z-axis. At a point (x, y, z) 
on the vortex core, the phases of the four surrounding 
points (x + Ax, y, z), [x, y + Ax, z), [x — Ax, y, z), and 
(x,y — Ax, z), all in the x — y plane, rotate by 2-k. By 
extending this analysis to the y — z and z — x planes, we 
can find all points on a vortex core. These two meth- 
ods of visualizing vortex cores are completely equivalent 
numerically, which means that quantized vortices are def- 
inite topological defects and thus different from eddies in 
classical fluid. An example of a fully developed turbu- 
lence with tangled vortices is shown in Fig[7| (a). The 
system at this time (t — 6) is in a quasi-steady turbu- 
lent state. Moreover, the energy spectrum e 2 / 3 k~ v is in 
quantitative agreement with the Kolmogorov law (Fig. [7| 
(b)). We thus conclude that the QT in quasi-steady state 



has the same statistics as the steady CT that follows the 
Kolmogorov law. The inertial range Ak < k < 2-7r/£ is 




(a) (b) 

FIG. 7: Simulated configuration of quantized vortices and 
the energy spectrum at t = 6. (a) The configuration, (b) The 
energy spectrum, which was obtained by making an ensem- 
ble average of 20 initial states. The solid line refers to the 
Kolmogorov law in eq. 

now 0.2 < k < 6.3, which is wider than that in the sim- 
ulations by Araki et al. an d Nore et al. [20|. The 
Kolmogorov constant is estimated to be C = 0.32 ± 0.1, 
and smaller than that in CT which is estimated to be 
1-4 < C < 1.8 [13. Araki et al. got the Kolmogorov 
constant C = 0.7 in their numerical simulation of the 
vortex-filament model and this is also smaller than that 
in CT [l9|. This small Kolmogorov constant may be, 
therefore, characteristic of QT. 



VI. STEADY TURBULENCE 

By introducing energy injection to the QT through a 
moving random potential, we model steady QT. The ad- 
vantage of steady QT over decaying QT are the following. 
First, steady QT gives a clearer correspondence with the 
Kolmogorov law; this is because the original statistics in 
the inertial range has been developed for steady turbu- 
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lence. Second, it enables us to confirm the presence of 
the energy-containing range, the inertial range and the 
energy-dissipative range of QT, which also exist in CT. 
Third, in all ranges, we can obtain the time-independent 
energy flux in the wave number space. Consequently, it 
is possible to reveal the cascade process of QT, which 
occur only by quantized vortices and short-wavelength 
excitations. This is later shown in Fig. [5J 

For energy injection at large scales, we introduce the 
moving random potential V(x, t), 

[i-7(sM)]J^*(!M) 

= [-V 2 - »(t) + g\<t>(x,t)\ 2 + V(x,t)]$(x,t). (26) 

We can numerically calculate the time development of 
the GP equation l|26[) by the same procedure that is dis- 
cussed in Sec. IIIII except h(k,t) in eq. (|T§|> is replaced 
by the Fourier transformation of h(x,t) = [g\$>(x, t)\ 2 + 
V(x,t)]$(x,i). To make the moving random potential 
V(x, t), we used a similar method as that used for making 
the random phase <f>o(x) in Sec. 0as the initial state of 
the wave function $>(x, t = 0)=fo exp[icpo(x)]. We place 
random numbers between and Vq in space-time (x,t) 
at every distance X for space and T for time, and con- 
nect them smoothly by using the four-dimensional spline 
interpolation [2^. We then obtain the moving random 
potential that has the following Gaussian two-point cor- 
relation 



(V(x,t)V(x',t')) 



Vq exp 



(x 



Y {t-t'f 



2X1 



2T* 



(27) 



This moving random potential has the characteristic spa- 
tial scale Xq and thus quantized vortices of radius Xq 
are nucleated when Vq is strong enough. We define the 
wave number separating the energy-containing range and 
the inertial range as 2tt/Xq. The wave number 27r/£ be- 
tween the inertial range and the energy-dissipative range 
is defined by the dissipation term y(x,t) the same as 
that in decaying turbulence. Therefore, our steady QT 
has the energy-containing range k < 2tt/Xq, the iner- 
tial range 2tt/Xq < k < 2ir/£ and the energy-dissipative 
range 2ir/£_ < k. 

We now show the result of our numerical simula- 
tion of steady turbulence for Vq = 50, Xq = 4, and 
To = 6.4 x 10~ 2 which are suitable numerical parame- 
ters to create fully developed turbulence. We start from 
the uniform wave function Q(x, t = 0) = fa = 1. Figure|H| 
shows the time development of the total energy E(t), the 
kinetic energy -Eki n (£), the compressible kinetic energy 
££ in (i), and the incompressible kinetic energy E^ in (t). 
Both figures show that the incompressible kinetic energy 
E^ in (t) is always dominant in the total kinetic energy 
E)dn(t)'i the moving random potential contributes to the 
nucleation of vortices rather than that of sound waves. 
As shown in Fig. [51(b), the four energies are almost con- 
stant for t > 25; hence, we can obtain a steady QT after 
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FIG. 8: Time development of E(t), E kin (t), E^ ln (t) and 
E k i n (t) at (a) the initial stage < t < 5 and (b) a later 
stage 25 < t < 30 (b). 



t ~ 25. The vortex configuration of steady QT looks 
similar to that of decaying QT in Fig. [7| (a). 

This steady turbulence is sustained by the competi- 
tion between energy injection and energy dissipation, 
which also results in a steady distribution of vortices. 
During steady QT, the energy flow in the wave num- 
ber space looks like that in Fig. The upper half of 
the diagram shows the kinetic energy E^ in (t) of vortices 
and the lower half shows the kinetic energy E^ in (t) of 
compressible excitations. In the energy-containing range 



energy-containing 
range 



2^ 

cucfgy-rlissipativti 
range 




Decay ol' 
vortices through 
recomicfu'oBft 



FIG. 9: The flow of the incompressible kinetic energy 
i?ii n (t) (upper half of diagram) and compressible kinetic en- 
ergy E^ in (t) (lower half) in wave number space. 

k < 2it/Xq, the system receives incompressible kinetic 
energy from the moving random potential with the nu- 
cleation of large quantized vortex loops. Throughout the 
Richardson cascade process of these quantized vortices, 
the energy flows from small to large k in the inertial 
range 2n/X < k < 2tt/£_. In the energy-dissipative 
range 2tt/£ < k, the incompressible kinetic energy trans- 
forms to compressible kinetic energy through reconnec- 
tions of vortices and the disappearance of small vortex 
loops. The moving random potential also creates long- 
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wavelength compressible sound waves that are another 
source of compressible kinetic energy and also produce 
an interaction with vortices. However, the sound waves 
have a very weak effect because E^. in (t) is much larger 
than £J£j n (t) as shown in Fig. [SJ moreover, sound waves of 
long wavelength have small scattering cross-section with 
vortex cores of the order of £. 

To further test the scenario in Fig. [§| we investigate 
the energy dissipation rate e of El in (t) and the flux of 
energy Il(fc, i) through the Richardson cascade in the in- 
ertial range. The energy dissipation rate e of E^lt) in 
a steady QT can be equated to e = —dEl in (t) /dt after 
switching off the moving random potential. This is be- 
cause the incompressible kinetic energy ££ in decays to 
the energy of compressible short-wavelength excitations. 
Figure^l (a) shows one example of the time development 
of E^. in (t) after switching off the moving random poten- 
tial. At the beginning, E^ in (t) decreases almost linearly 
at a rate from which we obtained e ~ 12.5 ± 2.3. On the 
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FIG. 10: (a) Time development of the incompressible kinetic 
energy E^ in (t) after switching off the moving random poten- 
tial in a steady QT. (b) Dependence of the energy flux II(fc, t) 
on the wave number k. The energy flux was calculated from 
an ensemble average of 50 randomly selected states at t > 25. 
The solid line is the energy dissipation rate e calculated from 
Figure (a). 



other hand, the energy flux H(k,t) of the incompressible 
kinetic energy in QT can be calculated by considering 
the scale-by-scale energy budget equation as well as CT 
0. To do this, we start from eq. H15[) with the moving 
random potential V(x,t): 

f(x,t)^(x,t) 

= j(x,t)[2Vf(x,t) ■ V0(x,t) + f(x,t)V 2 <j>(x,t)] 
+ V 2 /0M)-/0M)[V#r,t)] 2 
- [gf(x, tf - /i(t) + V(x, t)]f(x, t), (28) 

and again neglect the terms Vf(x, t) and df(x, t)/dt be- 
cause they are not effective at scales exceeding £. In- 
troducing a new variable p(x,t) = f(x,t)'V(p(x,t), this 



equation reduces to 




at 



p(x, t) + p(x, t) ■ Vv(x, t) 
= 7 {x,t)V[V-p(x,t)] + F(x,t), 



(29) 



where F(x,t) = -V[gf(x,t) 2 - + V(x,t)]f{x,t). 
Next, for any function A(x), we introduce the 
low-pass and high-pass filtered functions Af(x) = 

T, q < k Mq)eMi<l- x } an d A>(x) = J2 q>k Mq) ex P[i<? ' x ] 
respectively, and the operator of the low-pass filter L k : 
A{x) — > Af (x). Applying L k and taking the scalar prod- 
uct with p5(a:,f) to eq. (|29|l . we obtain 

l§- t \Pk (*> *)] 2 + Pk (*,*)• L k [p(x, t) ■ Vv(x, t)} 



1 (x,t)p<(x,t) -V[V-p< (x, t)} 
+ p<(x,t)-F<(x,t). 



(30) 



Extracting the incompressible part of eq. i|3U|) and inte- 
grating it over space yields the scale-by-scale incompress- 
ible kinetic energy budget equation 

J^(fc, t) + n(fc, t) = -Q(k, t) + &(k, t). (3i) 

Here we introduce the cumulative incompressible kinetic 
energy #(fc, t) between and k 



£{k,t) = - J dx [{p<(x,t)Y} 2 , 
the cumulative energy dissipation f2(fc,£) 



(32) 



0(fc,*) = -| J dxj(x 7 t){Pk(x,t) -V[V • P ^x,t)]}\ 

(33) 

the cumulative energy injection JP(k, t) 



2 

N 



dx {p<(x,t)-F<(x,t)Y, 



(34) 



and the energy flux through k of the incompressible ki- 
netic energy 



IL(k,t) 
2 

~ N 



dx[{p<(x,t)-L k \p(x,t) ■ Vv(x,t)]Y]. (35) 



Equation l|31|l can be interpreted as follows: at a given 
scale k, the rate of change of the incompressible kinetic 
energy is equal to the energy injected by the force &{k, t) 
minus the energy dissipation 17(/c,t) minus the energy 
flux n(fc,i) to smaller scales. The calculated energy flux 
n(fc, t) in Fig. E3(b) is nearly constant at H(k, t) ~ 11±1 
and is consistent with the energy dissipation rate e in 
the inertial range. These results indicate that the incom- 
pressible kinetic energy steadily flows in wave number 
space through the Richardson cascade at the constant 
rate n, and finally dissipates to compressible excitations 
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at the rate e ~ II. This energy flow is shown in the dia- 
gram of Fig. [5J The energy flux II is slightly smaller 
than the energy dissipation rate e, which may be at- 
tributable to the weak interaction between vortices and 
sound waves. 

To complete this picture of QT, we calculated the en- 
ergy spectrum E^ in (k). The result is plotted in Fig. ^2 
Because the case is similar to that of decaying turbulence, 
the energy spectrum is quantitatively consistent with the 
Kolmogorov law in the inertial range 2tt/Xq < k < 2ir/£, 
which is equivalent to 0.79 < k < 6.3. This inertial range 
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FIG. 11: Energy spectrum E^ in (k,t) for quantum turbu- 
lence. The plotted points are from an ensemble average of 
50 randomly selected states at t > 25. The solid line is the 
Kolmogorov law Ce 2/3 k~ 5/3 . 

is slightly narrower than that of the decaying turbulence 
because the moving random potential sets the energy- 
containing range. The resulting Kolmogorov constant C 
is C — 0.55 ± 0.07, which is smaller than that in CT 
as well as that for decaying turbulence. These smaller 
Kolmogorov constants C in both decaying and steady 
superfluid turbulence than that in CT may be charac- 
teristic of QT. We propose clarifying that as one of our 
future works. 



VII. CONCLUSION 

In this paper, we argue that QT can be a prototype of 
CT. To support this argument, we used a modified CP 



equation to investigate the dynamics of QT and com- 
pared the results to those of CT. In the original GP 
equation, turbulence is suppressed by compressible short- 
wavelength excitations. To overcome this problem, we 
introduced a dissipation term that has effects only at 
scales smaller than the healing length and then numeri- 
cally studied the modified GP equation. The main results 
are as follows. 

(i) The dissipation term does not disturb the dynamics 
of quantized vortices and enables one to investigate 
the statistics of quantized vortices in QT free from 
short- wavelength excitations. 

(ii) For decaying turbulence that starts from a random 
phase configuration of the wave function, the dissi- 
pation is effective for compressible excitations. This 
results in a incompressible kinetic energy spectrum 
that is consistent with the Kolmogorov law. 

(iii) A steady QT can be made by introducing energy 
injection as a moving random potential. In both 
steady and decaying turbulence, the incompressible 
kinetic energy dominates the total kinetic energy 
because of the dissipation term. The energy flux 
of the incompressible kinetic energy has a constant 
value in the inertial range and is consistent with the 
energy dissipation rate at small scales. The energy 
spectrum of the incompressible kinetic energy is also 
quantitatively consistent with the Kolmogorov law 
in the inertial range. These two results support the 
picture of QT shown in Fig. [§1 

These results indicate that QT is similar to CT and thus 
QT can become a prototype of turbulence for under- 
standing the inertial range, the Kolmogorov law, and the 
Richardson cascade process. We hope that these results 
will lead to new insights into the nature of turbulence. 
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